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We compute non-perturbatively the renormalization coefficients of scalar and pseudoscalar operators, local 
vector and axial currents, conserved vector and axial currents, and O^E^ 2 over a wide range of energy scales 
using a scaling technique that connects the results of simulations at different values of coupling p. We use the 
domain wall fermion formulation in the quenched approximation at a series of three values of /3 = 6.0, 6.45, 7.05 
corresponding to lattice spacing scaling by factors of two. 



1. INTRODUCTION 

In lattice QCD a method for full non- 
perturbative renormalization of operators by 
means of Monte-Carlo simulations was proposed 
m 0. To extend the range of momenta over 
which operators are renormalized, we combine 
simulations at different values of coupling con- 
stant (3 but with overlapping regions of phys- 
ical momenta, relating renormalization coeffi- 
cients between momenta scales from 2 to 10 GeV. 

2. NONPERTURBATIVE RENORMAL- 
IZATION METHOD 

The renormalized operator 0(/i) is defined as 



O(ju) = Z (p;a)O bare (a) 



(1) 



where for an operator with n quark fields Zo (p] a) 
is defined by the renormalization condition 



ZoifJ-; a)Z^ n/2 (^ a)T {p; a)\ p 2 



1 



(2) 



Here To {p) is a truncated Greens function traced 
with a projector Po on the tree-level opera- 
tor. It is calculated with far off-shell exter- 
nal momenta in Landau gauge. For quark bi- 

linears, To(p;a) — Tr ( h-o{p\ a )Po ) /12, where 
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A (p;a) = S(p;a) 1 G (p;a)S(p;a) \ and 
Go (p', a) is non-amputated Greens function. 

The quark field renormalization coefficient is 
defined through the Ward identity for the con- 
served vector current, 



Z ^,a) = ^Tr ( 7 "A£ c (p; a )) p2= 



(3) 



3. SCALING TECHNIQUE 



When scaling violations are small, the renor- 
malization coefficients calculated at two differ- 
ent values of the lattice spacing are related to 
each other by an overall multiplicative coefficient 
which is independent of the physical conditions 
such as physical volume V or momentum scale //, 



Zp(p, Vjaf) 
Z (v,V;a) 



= R(a',a) 



(4) 



For each operator we perform a series of simu- 
lations with lattice spacing a increasing by factors 
of two. The number of lattice sites in all direc- 
tions is fixed. The physical volume is finite and 
different in each simulation, increasing by factors 
of 2 4 . 

We rescale the Z-factors using property (H) so 
that they are all defined at a single value of (3 — 
7.05. To compute the scaling coefficients R(a', a), 
a series of additional simulations is used with half 
the number of the lattice sites in each direction 
and twice as large lattice spacings, so that the 



2 



physical volume and momenta match those in the 
main series. We use fits quadratic in fia to remove 
the discretization error corrections to (Q) that for 
domain wall fermions are of order (fia) 2 , 

Z {ii,V;a) 



R(a,2a) = 



0((pa) 2 



(5) 



Z {^V;2a) 

Once we know the rescaling coefficients for pairs 
of lattice spacings, R(a,2a), R(2a,4a), etc., we 
use them to rescale the Z-factors obtained at dif- 
ferent /3's, 

Z (n,V;a)=R(a,2a)Z (n,V;2a) (6) 
Z {li, V; a) = R{a, 2a)R(2a, 4a)Z (fJ>, V; 4a) (7) 
etc. 

An important advantage of this method is that 
the physical volume is explicitly finite in the sim- 
ulations at intermediate /3's. The infinite volume 
limit can be taken at the largest (3 to match to 
perturbation theory. 
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Figure 1. Z^/Zy in the chiral limit for the local 
vector current in the same physical volume and 
rescaling coefficient i?(a(6.45), a(7.05)). 

4. NUMERICAL RESULTS 

In our simulations we use lattice volumes 8 4 
and f6 4 at three values of/3 = 6.0, 6.45, and 7.05 
with the inverse lattice spacing of 1.96, 3.92, and 
7.84 GeV respectively §]. 
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Figure 2. Z^/Z 
7.05. 
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at different /3's, rescaled to (3 



4.1. Local Vector and Axial Currents 

Figure [l] illustrates the procedure of obtaining 
the rescaling coefficient for the local vector cur- 
rent from two simulations at different /3's in the 
same physical volume. The two sets of symbols 
on the main graph correspond to Z^/ZvifJ-, V; a) 
for two values of the lattice spacing at (3 = 7.05 
and 6.45 different by a factor of two. The smaller 
plot shows a linear extrapolation of the ratio of 
these Z-factors to [i = according to (||) . 

Figures || and || show rescaled Z-factors defined 
at /3 = 7.05. 
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Figure 4. Z^/Zs at different /3's, rescaled to (3 = 
7.05. Solid line - RG running of Z^jZ s . 

4.2. Scalar and Pseudoscalar Densities 

When taking the chiral limit of scalar and pseu- 
doscalar densities at energy scales of a few GeV 
in finite volume, one needs to remove non-leading 
1 jp 2 contributions which have a singular behavior 
as m — ► 0. These 1/p 2 contributions are propor- 
tional to the chiral condensate, 



r>b 2 ) = fV) + ^ 

Zp p z rrif 



z ± 
z s 



(p 2 ) 



C d{4>ip) 
v 2 drrif 



(8) 



(9) 



For very small masses the zero modes make a 
significant contribution to {ipip) with a term in- 
versely proportional to m/. This term con- 
tributes to @ as —C'/{m 2 p 2 ). After subtracting 
these terms, the scaling procedure is applied with 
the results shown in Figures | and |. 

4.3. Conserved Vector and Axial Currents 

The conserved currents can be used to compute 
Z^ from (|). With the DWF action, 

E - gi E i 101 

s=L a /2 s=0 



1 — 7 

jp(x,s) = ip s (x) — ^-Ufj,(x)ip s (x + fx) 
-$ a ( x + ii)±^Ut(x)il> a (x), 



(11) 
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Figure 5. Z^/Zp at different /3's, rescaled to (3 = 
7.05. Solid line - RG running of Z^/Zp. 

where <jy = +1 and a a = — 1. 

Due to the nonlocal character and summation 
over s, these matrix elements can be very expen- 
sive to calculate. To reduce the amount of com- 
putational time, we used a random source esti- 
mator to compute the part of the sum between 
s = 1 and L s — 2 and used only one component 
of (p) for different momenta to calculate 

Tr hT^p)) = Tr{ 7 °r°(p) + 7 r (p° p 1 ) 

-°rV «-» p 2 ) + 7 °rV ~ p 3 )} . (12) 



+7 U 
We set m/ 



0.05 in these calculations. Figures 
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Figure 6. Z,p from conserved vector current, 
rescaled to (3 = 7.05. Solid line - RG running. 
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Figure 7. Z^ from conserved axial current, 
rescaled to (3 = 7.05. Solid line - RG running. 

|^ and []] show rescaled defined at j3 = 7.05 ob- 
tained from the conserved vector and axial cur- 
rents. 

4.4. B K 

Bk is a phenomenologically important param- 
eter relevant to the calculation of K°—K° mixing, 



Bk(jj) 



{K°\O£t 2 M\K ) 

8 f 2 2 



(13) 



The renormalization of Bk is the same as that of 
Oll =2 = (^7l^) 2 - Using the upper indices for 
color and the lower for spin, the non-amputated 
Fourier transformed four-point Greens function 
with equal off-shell external momenta can be 
written as 



AR ^ hsS(p\o)H 



(14) 

where r^(p)^ = S(p\Q)i?H„ p MCpWtsC • 
S(p\0) is a propagator on a single configuration. 
The projection operation is defined as 



(p; a) 



V V Ad (n\ AABB 



(15) 



32N C (N C + 1) 
with the amputated Greens function given by 
A BK (P)^ = S-\p) AA 's~\p)CC' 

■ GBAP^sP'S-'iP^ffS- 1 ^ . (16) 
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Figure 8. ^ 
f3 = 7.05. Solid line 



Z?i./Z as=2 at different /3's, rescaled to 



RG running. 



Applying our earlier analysis then gives the re- 
sults shown in Figure ||. 

We conclude that this combination of recursive 
scaling and NPR offers a promising connection 
between the perturbative and non-perturbative 
regimes-a connection that is required if pertur- 
bative errors are to be properly controlled. The 
results reported here suggest a consistent picture 
for the scaling of renormalization factors over a 
momentum range between 2 and 10 GeV for the 
case of quenched domain wall fermions. More ex- 
tensive calculations will be required for a com- 
plete demonstration of this approach. 
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